Shortcut to adiabaticity for an interacting Bose-Einstein condensate 
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We present an investigation of the fast decompression of a three-dimensional (3D) Bose-Einstein 
condensate (BEC) at finite temperature using an engineered trajectory for the harmonic trap- 
ping potential. Taking advantage of the scaling invariance properties of the time-dependent Gross- 
Pit aevskii equation, we exhibit a solution yielding a final state identical to that obtained through 
a perfectly adiabatic transformation, in a much shorter time. Experimentally, we perform a large 
trap decompression and displacement within a time comparable to the final radial trapping period. 
By simultaneously monitoring the BEC and the non-condensed fraction, we demonstrate that our 
specific trap trajectory is valid both for a quantum interacting many-body system and a classical 
ensemble of non-interacting particles. 

PACS numbers: 37.10.-x, 67.85.-d 



Quantum adiabatic transformations [H, 0, in which 
the system's parameters are varied slowly enough such 
that no transition between instantaneous eigenstates oc- 
cur, play a central role in physics. For instance, schemes 
based on adiabatic passage have been proposed to pre- 
pare non classical states 0, Q or to produce new strongly 
correlated states [5]. Quantum adiabatic computation 
is attracting a lot of attention 0, Q. Adiabatic trans- 
formations using various time-dependent potentials are 
routinely performed in experiments on ultracold gases. 
However adiabatic transformations are typically slow [i| , 
while practical applications and experimental constraints 
such as finite lifetime or coherence time [9] require faster 
processes. 

This contradiction motivated the search for rapid 
schemes reproducing or approaching ideal adiabatic tran- 
sitions. "Exact" methods [Io|, El, here referred to 
as shortcuts to adiabaticity, yield a final state strictly 
identical to that obtained via an adiabatic transforma- 
tion, while other approaches [l2l-[l5| use minimization 
techniques to optimize the transition to a target state. 
Among the former, some strategies require introducing 
additional terms in the Hamiltonian [To[, Efij j others con- 
sist in engineering the time dependence of the parameters 
to avoid unwanted transitions [17]. In spite of the large 
theoretical literature, few experiments were conducted on 
classical systems 18-2o|, and even fewer in the quantum 
regime [21] . 

In this Letter, we perform the rapid shortcut de- 
compression of a 3D interacting BEC confined in an 
anisotropic harmonic trap. The trap frequencies are de- 
creased by a factor of 9 (radially) and 3 (axially) in a 
time comparable to the final radial trapping period, us- 
ing a trajectory based on the scaling properties of the 
time-dependent Gross-Pitaevskii equation (GPE) in the 
Thomas- Fermi (TF) limit [22]. This shortcut trajectory 
leads to a final state identical (in theory) to the equi- 
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FIG. 1: Linear versus shortcut BEC decompression. We com- 
pare the time evolution of the BEC after two different decom- 
pression schemes: (A) a 30-ms-long linear ramp and (B) the 
shortcut trajectory (see text). The center of mass motion has 
been subtracted from these time of flight images for clarity. 



librium state obtained via a perfectly adiabatic process. 
Experimentally, we demonstrate that the collective exci- 
tations [H, Q associated with the rapid trap decompres- 
sion are strongly reduced by our shortcut scheme (Fig. [I]), 
the residual excitation being due to experimental imper- 
fections. Furthermore, we show that the trajectory is 
also valid for a classical ensemble of non-interacting par- 
ticles, as demonstrated by monitoring the non-condensed 
fraction of the finite-temperature BEC. 

We describe our system by a zero-temperature BEC 
plus a thermal cloud, assumed to behave independently. 
The BEC component thus obeys the 3D GPE 
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where -0(r, £) is the wave function of the condensate, m 
the mass, N the number of particles, and U = 47rh 2 a s /m 
the interaction coupling constant (a s is the s-wave scat- 
tering length). The time-dependent trapping potential 
has a cylindrical symmetry along the horizontal axis (y) 

U(r,t) = ^mu ± (t) 2 (x 2 +z 2 ) + ^mu 2 (t)y 2 + mgz, (2) 

with initial and final angular frequencies cj {±,||} and 
= ^o{J_, || > / ||' respectively. It is worth stress- 
ing here that decreasing the trap frequencies not only de- 
compresses the BEC but also translates the harmonic po- 
tential minimum vertically by Az = — g(l/ujQ ± — 1/uo 2 ^). 
The objective is to engineer a trajectory (t) connect- 
ing the equilibrium states in the initial and final poten- 
tials. We stress that the BEC is not at equilibrium at 
any time during the trajectory, but only at t = and 

t = t f . 

Equation (pQ) is invariant under the scaling and trans- 
lational transformation 

^(r,t) = (6i6||)- 1 / 2 X (p,r(t))exp[i0(r,t)] (3) 
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with p x = x/b ± , p y = y/b\\, p z = z/b± + ga/u$ ± and 
r(t) = Jq dt' /\b\(t')b\\ (£')], in the following situations: 

or (ii) for 



either (i) in the non-interacting limit [22, 
a suitable driving of the interaction term U via Feshbach 
resonances [25[; or (hi) in the Thomas- Fermi (TF) limit 
(i.e. neglecting the term V 2 x) [26[. 

In the latter case, one can derive three coupled differen- 
tial equations satisfied by the scaling (bj_i|) and shifting 
(a) parameters, the condition of initial and final equilib- 
rium imposing sixteen independent boundary conditions. 
Our procedure to engineer (t) is to reduce the dimen- 
sionality of the problem by looking for trajectories with a 
constant by (and thus a constant axial size of the BEC), 
yielding a decompression with 7^ = 7y • In this case, 6_l,|| 
and a should satisfy 

b±(t) + b ± (t)ujl(t) = wg ± /6 ± (t) 3 (4) 
u\\(t) =uj oll /b ± (t) (5) 

b±(t) 4 d(t) + 2b±(t) 3 b ± (t)a(t) + u>l ± a(t) - uj% ± b ± (tf = 0, 

(6) 

for which we can exploit the procedure outlined in [2o| . 
Since Eqs. (jlj and (J6j) are identical to those obtained for a 
thermal cloud [20], the shortcut trajectory exhibited here 
is valid both for the BEC (in all directions) and for the 
thermal fraction in the radial directions. Indeed, since 
Eq. dSJ) does not yield a fixed axial size for a thermal 
cloud, this shortcut trajectory will not work for the axial 
direction in the thermal case. 

Our ultracold 87 Rb atoms are trapped in a quadrupole- 
Ioffe-configuration (QUIC) magnetic trap [27j. In its 
compressed initial state this trap is anisotropic with 
radial and axial frequencies ujo±/2tt = 235.8 Hz and 
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FIG. 2: Shortcut BEC decompression in 30 ms. We plot 
the shortcut trajectories uj±{t)/2n (solid) and uj\\{t)/2i: (x5, 
dashed). The insert compares the subsequent evolution of the 
BEC's center of mass (open symbols) and radial size (solid 
symbols) for the shortcut (stars) and linear (squares) decom- 
pressions (GPE simulation). 



co>o||/27r = 22.2 Hz respectively. We can tune these inde- 
pendently by adjusting the QUIC current (affects both 
uj_\_ and cj||) and the current running through a pair of 
Helmoltz coils aligned along y (affects only ou±) [28|. In 
the present paper, we perform a 30-ms-long radial decom- 
pression of the trap by a factor of 9, yielding a final radial 
frequency ujf±/27T = 26.2 Hz. Because of the condition 
7_l = 7y imposed above, the axial frequency is reduced 
by a factor of 3 to a final value U0f\\/2ii = 7.4 Hz. If 
performed adiabatically, this decompression reduces the 
chemical potential by a factor of 9. The corresponding 
shortcut frequency trajectories are represented on Fig. [2j 
The insert shows the expected evolution of the vertical 
center of mass position of the BEC and its radial size 
as obtained from GPE numerical simulations. For the 
shortcut trajectory (stars), these quantities remain sta- 
tionary after the end of the decompression indicating that 
an equilibrium final state is reached. For comparison, we 
plot (squares) the same quantities for a linear ramp of 
same duration, yielding large dipolar and breathing os- 
cillations [29|. 

The experimental procedure is as follows. We first pro- 
duce a BEC by RF evaporation in the compressed trap. 
In the experimental runs presented here, the condensed 
fraction (N = 1.3 x 10 5 ) represents 60% of the total num- 
ber of atoms. The initial temperature, inferred from the 
size of the non-condensed fraction after time of flight, is 
To = 130 nK. We then apply a decompression sequence, 
hold the ultracold cloud for a certain time th in the de- 
compressed trap, then release it and monitor the cloud's 
parameters after a 28 ms time of flight via absorption 
imaging. This time of flight is close to the critical time 
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30 ms) where the aspect ratio of the decompressed 
BEC inverts, which explains its isotropic aspect on Fig.[TJ 
By varying t^, we characterize the magnitude of the var- 
ious modes excited by the decompression process. To 
extract quantitative estimates, we fit the 2D column den- 
sity profiles by a two-component distribution allowing for 
different angles (see discussion below) for the BEC (TF 
profile) and thermal fraction (Gaussian profile). The fit 
results are averaged over three different images taken in 
the same conditions. Throughout this paper we will com- 
pare three different decompression schemes: an abrupt 
jump from the initial to final frequencies, a linear ramp 
of duration 30 ms, and the shortcut trajectory depicted 
on Fig. El 

Fig. [U illustrates the efficiency of our shortcut method. 
The vertical axis on the figure corresponds to the direc- 
tion of gravity (z), and the horizontal one to the axial 
direction (y). The field of view of each image is 576 jam 
x 576 /im, the indicated value corresponding to the time 
th spent in the decompressed trap. The center of mass 
motion has been subtracted from these images for clarity 
(see Fig. 03(A)). (A) corresponds to a 30-ms-long linear 
decompression while (B) is obtained with the shortcut 
trajectory. Qualitatively, we observe that in the first 
case the BEC undergoes large deformations character- 
ized by an oscillatory behavior of its aspect ratio. More 
unexpectedly, the BEC is also seen to oscillate angularly 
about the horizontal axis, reflecting the excitation of a 
scissors mode [30, 

El- 

This "parasitic" excitation will be 
discussed at the end of the paper. Applying the shortcut 
scheme results in a drastic suppression of the BEC de- 
formations (Fig. QjB)). Indeed, the BEC remains almost 
stationary throughout the whole 170 ms range after de- 
compression, very close to the targeted equilibrium state. 

Fig. [31(A) shows the oscillation occurring at ujf± of 
the cloud's center of mass position along the vertical (af- 
ter a 28 ms time of flight). This dipole mode is excited 
by the large vertical displacement (357 /im) of the trap 
center. The circles, squares and stars correspond to the 
abrupt, linear and shortcut decompressions respectively. 
The lines are sinusoidal fits. The shortcut decompression 
reduces the amplitude of the dipole mode by a factor of 6 
and 4.3 when compared to the abrupt and linear schemes 
respectively. The corresponding peak vertical velocity is 
1.8 recoil. The origin of this residual dipole oscillation 
will be discussed at the end of the paper. Note that the 
center of mass oscillation is totally uncoupled from the 
internal degrees of freedom of the BEC, even in the pres- 
ence of interactions [291 ] . 

We report on Fig.[3](B) the evolution of the aspect ratio 
of the BEC (again after a 28 ms time of flight). For the 
abrupt and linear schemes, it presents strong oscillations 
at a frequency of 47 Hz which is consistent with a radial 
breathing mode of frequency « 2ujfj_/27r 0, 0, [HI . A 
Fourier transform analysis of the BEC size oscillations 
also reveals the presence of an axial breathing oscilla- 
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FIG. 3: (Color online) Decompression-induced excitations. 
We report the temporal evolution of (A) the center of mass 
position and (B) the aspect ratio of the BEC after three differ- 
ent decompression schemes: an abrupt decompression (open 
circles); a 30 ms linear ramp (squares); the 30ms shortcut 
trajectory (stars). All measurement are performed after a 28 
ms-long time of flight. 



tion at 12.5 Hz (« ^Jhj2 u f \\/2it |29|). The specific 
shape of the aspect ratio fluctuations on Fig. [3(B) re- 
sults from the interplay between the various modes and 
depends on the decompression trajectory as confirmed by 
GPE simulations. Using the shortcut trajectory (stars on 
Fig. 119(B)) strongly inhibits these breathing-like excita- 
tions: the standard deviation of the aspect ratio varia- 
tions versus th is reduced by a factor of 12 and 10 when 
compared to the abrupt and linear decompressions re- 
spectively. After the shortcut decompression and the 28- 
ms-long expansion, the BEC has an average TF radius 
of 46.8 jam (close to the GPE prediction of 43 jam) with 
a standard deviation below 12% (see Fig. IU(A)). 

Fig. [4] illustrates the impact of the shortcut decom- 
pression on the thermal fraction. Here, as mentioned be- 
fore, the shortcut trajectory is expected to work only for 
the radial directions. Indeed, we observe on Fig. IU(B) a 
breathing oscillation at 2 x ctjy || / 2tt of the thermal cloud 
along the axial direction (open circles on lower panel), 
while its size along the vertical (filled circles) remains 
stationary, corresponding to a temperature Tf = 22 nK 
(a factor of 5.9 below the initial temperature). In the 
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FIG. 4: BEC versus thermal cloud decompression. We plot 
the sizes of the BEC (A) and thermal component (B) versus 
th for the shortcut trajectory. The filled and empty sym- 
bols correspond to the radial (vertical) and axial directions 
respectively. 



case of the linear ramp, we observe breathing oscillations 
along both axial and radial dimensions. For compari- 
son, we show on (A) the sizes of the BEC along both 
axis, which are stationary since the shortcut trajectory 
is valid both radially and axially. 

We now discuss the various imperfections of the ex- 
periment, which limit the performances of our shortcut 
decompression and result in a final state different from 
the targeted equilibrium one. A first problem is the mis- 
match between the theoretical and experimental trap fre- 
quency trajectories. For instance, we measure from the 
dipole oscillation Uf±/27r = 23.6 Hz which is 10% below 
the targeted value of 26.2 Hz. Since the cloud's dynam- 
ics is quite sensitive to the final stage of the decompres- 
sion, this can account for the residual dipole oscillation of 
Fig. 02 A). Another possible issue is the trap anharmonic- 
ity, which is enhanced in the decompressed trap due to 
gravity. The impact of the anharmonicity could be re- 
duced by using a smoother trajectory [18|,[20|. The trap's 
non-ideal geometry is also responsible for the excitation 
of the scissors mode observed on Fig. HJ In Ref. [3l|, an 
angular velocity was imparted to the BEC by suddenly 
tilting the trap, and the subsequent monochromatic an- 

Jm was used as an evidence 



gular oscillation at 

for superfluidity. Here, the angular momentum is com- 
municated to the cloud during the trap decompression 



because the trap eigenaxes tilt slightly as the trap cen- 
ter moves downwards due to gravity. This results in an 
angle of 3° between the axial directions of the initial and 
final traps in the vertical plane. Note that the large am- 
plitude of the scissors oscillation on Fig. Q] is due to the 
magnification effect of the time of flight [33|, [34| . In-situ 
measurement show an amplitude of the scissors oscilla- 
tion compatible with the trap tilt angle of 3°. 

In conclusion, we presented in this paper a method 
to perform shortcut-to-adiabaticity transformations on a 
3D interacting BEC, using a specifically designed param- 
eter trajectory for the harmonic trapping potential. The 
performances could be further improved using better- 
controlled potentials such as in optical traps or lattices, 
where time-dependent manipulations are also easier and 
faster. Very short transition times could in principle be 
achieved by transiently applying negative (i.e. expelling) 
curvatures Further work may include the direct 

comparison with other methods such as "bang-bang" [35| 
or optimal control techniques. More general shortcut so- 
lutions will also be searched for, and applied to other 
dimensionalities or non-harmonic potentials [37j. These 
fast transition methods are not restricted to cold atom 
manipulation, and can be readily adapted to topics as 
diverse as e.g. macroscopic resonator cooling [36|, tem- 
poral [l6| and spatial [38[ coherent population transfer, 
or quantum computation [7]. 
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